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APPARATUS FOR GENERATING FOCUSED 
ELECTROMAGNETIC RADIATION 

The present invention relates to the generation of electromagnetic ra- 
diation and, more particularly, to an apparatus and method of generating 
focused pulses of electromagnetic radiation over a wide range of frequencies. 
More particularly it relates to an apparatus and method for generating pulses 
of non-spherically decaying electromagnetic radiation. 

The present apparatus and method axe based on the emission of electro- 
magnetic radiation by rapidly varying polarisation or magnetisation cxnrent 
distributions rather than by conduction or convection electric currents. Such 
cvirrents can have distribution patterns that move with arbitrary speeds (in- 
cluding speeds exceeding the speed of Ught in vacuo), and so can radiate over 
a much wider range of frequencies than their conventional counterparts. 

Furthermore, intensities of normal emissions decay at a rate of i2~^, 
where R is the distance from the source. It has been noted, however, that the 
intensities of certain pulses of electromagnetic radiation can decay spatially 
at a lower rate than that predicted by this inverse square law (see Myers et o/., 
Phys. World, Nov. 1990, p. 39). The new solution of Maxwell's equations set 
out below, for example, predicts that the electromagnetic radiation emitted 
from superluminally, circularly moving charged patterns decays at a rate 
of R"^. Another example is the electromagnetic radiation emitted from 
superlimiinally, rectilinearly moving charged patterns which decays at a rate 
of Ji-t. 

It will be widely appreciated that being able to employ such emissions 
for signal transmission, amongst other applications, would have significeuit 
commercial value, given that it would enable the employment of lower power 
transmitters and/or larger transmission ranges. The present invention pro- 
vides a method and apparatus for generating such emissions. 

According to the present invention there is provided an apparatus for 
generating electromagnetic radiation comprising; 

a polarizable or magnetizable medium; and 

means of generating, in a controlled manner, a polso-isation or mag- 
netisation ciurent with a rapidly moving distribution pattern such that the 
moving source in question generates electromagnetic radiation. 

The speed of the moving distribution pattern may be superluminal 
so that the apparatus generates non-spherically decaying electromagnetic 
radiation. 

The apparatus may comprise a dielectric substrate, a plurahty of elec- 
trodes positioned adjacent to the substrate, and the means for applying a 
voltage to the electrodes sequentially at a rate sufficient to induce a polarised 



region in the substrate which moves along the substrate with a speed either 
approaching or exceeding the speed of light. The dielectric substrate may 
have either a rectilinear or a circulsa- shape. 

The wavelength of the generated electromagnetic radiation may be in 
any range from the radio to a minimum deternained only by the limit to the 
scale of the structure in the dielectric medium (potentially optical, ultraviolet 
or even x-ray). 

Examples of the present invention will now be described with reference 
to the accompanying drawings, in which: 

Figure 1 is a diagram showing the wavefronts of the electromagnetic 
emission from a particular volume element (source point) S within the mov- 
ing polarised region of the polarizable medium of the present invention; 

Figure 2 is a graph showing the vsdue of a function representing the 
emission time versus the retarded position for differing somrce points a, b, c 
within the polarizable medirnn in question; 

Figinre 3 is a perspective view of the envelope of the wavefronts shown 
in Fig. 1; 

Figure 4 is a view of the cusp curve of the envelope shown in Fig. 3; 

Figm*e 5 is the locus of the possible source points which approach the 
observation point P along the radiation direction with the wave speed at the 
retarded time, a locus that is henceforth referred to as the bifurcation smrface 
of the observer at P; 

Figure 6 is a view of the cross sections of the bif\u:cation surface and 
the source distribution with a cylinder whose axis coincides with the rotation 
axis of the source; and 

Figures 7(a) and 7(b) are views of two examples of the apparatus of 
the present invention showing the dielectric substrate, the electrodes and a 
superluminally moving polarised region of , the dielectric substrate. 

Prior to description of the invention, it is appropriate to discuss the 
principles imderlying it. 

Bolotovskii and Ginzbm-g (Soviet Phys. Usp. 15, 184, 1972) and Bolo- 
tovskii and Bykov (Sovet Phys. Usp. 33, 477, 1990) have shown that the 
coordinated motion of aggregates of charged particles can give rise to ex- 
tended electric charges and cinrrents whose distribution patterns propagate 
with a phase speed exceeding the speed of light in vacuo and that, once 
created, such propagating charged patterns act as sources of the electromag- 
netic fields in precisely the same way as any other moving soxu"ces of these 
fields. That these sources travel faster than light is not, of course, in any 
way incompatible with the requirements of special relativity. The superlu- 
minally moving pattern is created by the coordinated motion of aggregates 
of subluniinally moving particles. 



We have solved Maxwell's equations for the electromagnetic field that 
is generated by an extended source of this type in the case where the charged 
pattern rotates about a fixed axis with a constant angular frequency. 

There are solutions of the homogeneous wave equation referred to, inter 
a/ia, as non-diffracting radiation beams, focus wave modes or electromagnetic 
missiles, which describe signals that propagate through space with unexpect- 
edly slow rates of decay or spreading. The potential practic£tl significance of 
such signals is clearly enormous. The search for physically realizable sources 
of them, however, has so far remained unsuccessful. Our calculation pin- 
points a concrete example of the sources that are currently looked for in this 
field by establishing a physically tenable inhomogeneous solution of Maxwell's 
equations with the same characteristics. 

Investigation of the present emission process was originally motivated 
by the observational data on pulsars. The radiation received from these 
celestial sources of radio waves consists of highly coherent pulses (with as 
high a brightness temperature as 10^^ °K) which recur periodically (with 
stable periods of the order of 1 sec). The intense magnetic field (~ 10^^ G) 
of the central neutron star in a pulsar affects a coupling between the rotation 
of this St ax and that of the distribution pattern of the plasma surrounding it, 
so that the magnetospheric charges and currents in these objects are of the 
same type as those described above. The effect responsible for the extreme 
degree of coherence of the observed emission from pulsars, therefore, may 
well be the violation of the inverse square law that is here predicted by our 
calculation. 

The present analysis is relevant also to the mathematically similar prob- 
lem of the generation of acoustic radiation by supersonic propellers and he- 
licopter rotors, although this is not discussed in detail here. 

We begin by considering the waves that are emitted by an element of 
the superluminally rotating source from the standpoint of geometrical optics. 
Next, we calculate the amplitudes of these waves, i.e. the Green's function 
for the problem, from the retarded potential. We then specify the bifurcation 
sturface of the observer and proceed to calculate the electromanetic radiation 
arising from a superluminally moving extended source. The singuletrities of 
the integrands of the radiation integrals that occur on the bifm-cation surface 
are here handled by means of the theory of generalised functions: the electric 
and magnetic fields are given by the Hadamard's finite parts of the divergent 
integrals that result from differentiating the retarded potential under the 
integral sign. The theory is then concluded with a descriptive account of 
the analysed emission process in more physical terms and the description of 
examples of the apparatus. 



I. ENVELOPE OF THE WAVEFRONTS AND ITS CUSP 



Consider a point source (an element of the propagating distribution 
pattern of a volume sovirce) which moves on a circle of radius r with the 
constant angular velocity luB;,, i.e. whose path x(i) is given, in terms of the 
cylindrical polar coordinates (r,<^,z), by 

r = const., z — const., (p = (p + ijt, (1) 

where is the basis vector associated with z, and (p the initial value of (p. 

The wavefronts that are emitted by this point source in an empty and 
unbounded space are described by 

|xp-x(t)| = c(tp-t), (2) 

where the constant c denotes the wave speed, and the coordinates (xp,ip) = 
{rp,(pp,zp,tp) mark the spacetime of observation points. The distance R 
between the observation point xp and a source point x is given by 

|xp - x| = R{ip) = [izp - zf + rp^ + - 2rprcos(v?p -tp)]^, (3) 

so that inserting (1) in (2) we obtain 

R{t) = [(zp - z)^ + rp2 + - 2rprcos{(pp - (p - tjt)]^ = c{tp - t). (4) 

These wavefronts are expanding spheres of radii c{tp - 1) whose fixed centres 
(rp = Tytpp = (p + tjt, zp = z) depend on their emission times t (see Fig. 1). 

Introducing the natural length scale of the problem, cju), and using 
t = (^ - <p) /uj to eliminate t in favour of y?, we cein express (4) in terms of 
dimensionless variables as 

g = iP'-ipp + R{(p) = <f>, (5) 

in which R = Ruj/c, and 

<f>=(p — <pp (6) 

stands for the diflFerence between the positions (p = — ut of the source 
point and fpp = (pp - ujtp of the observation point in the (r, ^, 2:)-space. 
The Lagrangian coordinate (p in (5) Ues within an interval of length 27r (e.g. 
— TT < <^ < tt), while the angle (p, which denotes the azimuthal position of the 
source point at the retarded time t, ranges over (— oo, oo). 

Figure 1 depicts the wavefronts described by (5) for fixed values of 
{r,(p,z) and of 0 (or tp), and a discrete set of values of (p (or t). [In this 



figure, the heavier curves show the cross section of the envelope with the plane 
of the orbit of the source. The larger of the two dotted circles designates the 
orbit (at r = 3c/u;) and the smaller the hght cylinder (rp = c/u;).] 

These wavefronts possess an envelope because when r > c/ujy and so 
the speed of the source exceeds the wave speed, several wavefronts with 
differing emission times can pass through a single observation point simul- 
taneously. Or stated mathematically, for certain values of the coordinates 
(rp,(pp,zp;r,z) the function g{(p) shown in Fig. 2 is oscillatory and so can 
equal (f> at more than one value of the retarded position (p: a horizontal line 
(f> = constant intersects the curve (a) in Fig. 2 at either one or three points. 
[Fig. 2 is drawn for = 0, fp = 3, f = 2 and (a) z = zp, inside the 
envelope, (b) z = ic? on the cusp curve of the envelope, (c) z = 2zc — zp^ 
outside the envelope. The marked adjacent t\u-ning points of curve (a) have 
the coordinates {(p±, and ^©ut represents the solution of g{(p) = <f>o for 
a ^0 that tends to <f>^ from below.] 

Wavefronts become tangent to one another and so form an envelope 
at those points (rp,^p,zp) for which two roots of g{(p) — <f> coincide. The 
equation describing this envelope can therefore be obtained by eliminating (p 
between g — </> and dg/d(p = 0. 

Thus, the values of ip on the envelope of the wavefronts are given by 

dg/d(p =l-frp sm(¥?p - <p)/R{v>) = 0- (7) 

When the curve representing g{(p) is as in Fig. 2(a) (i.e. f > 1 and A > 0), 
this equation has the doubly infinite set of solutions <p ^ (p± + 2n7r, where 

= + 27r - arccos[(l T )/(ffp)], (8) 

A = (f|>-l)(f2-l)-(z-ip)^ (9) 

n is an integer, and (f,£;fp,£p) stand for the dimensionless coordinates 
rijj/c, zu/c, rpij/c and zpu/c, respectively. The function g{(p) is locally 
maximum at ^+ + 2n7r and minimum at (p- + 2n7r. 

Inserting (p = ip± in (5) and solving the resulting equation for 0 as a 
function of (rp, ip), we find that the envelope of the wavefronts is composed 
of two sheets: 

^ = 0± = 9(<P±) = 27r - arccos[(l =F )/(ffp)] + R±, (10) 

in which _ . 

= [{z - zpf + f2 + f|. - 2(1 T A^)]^ (11) 

are the values oi R ip = ip±. For a fixed source point (r, ^,z), equa- 
tion (10) describes a tube-like spiralling surface in the (rp, ^p, zp)-space of 



observation points that extends from the speed-of-light cylinder fp = 1 to 
infinity. [A three-dimensional view of the light cylinder and the envelope of 
the wavefronts for the same soiu*ce point (S) as that in Fig. 1 is presented 
in Fig. 3 (only those parts of these surfaces are shown which lie within the 
cylindrical volume fp < 9, —2.25 < zp — i < 2.25).] 

The two sheets <f> — (t>± of this envelope meet at a cusp. The cusp 
occiu's along the cxuA^e 

0 = 27r - arccos[l/(ffp)] + (f%f^ - 1)^ = (pc, (12a) 

i=ip±(f2,-i)§(f2-l)i (126) 

shown in Fig. 4 and constitutes the locus of points at which three different 
wavefronts intersect tangentially. [Figure 4 depicts the segment — 15 < fp — 
f < 15 of the cusp curve of the envelope shown in Fig. 3. This curve touches — 
and is tangent to — the light cylinder at the point (fp = l^zp = z^(f> = 
^c|ri»=i) on the plane of the orbit.] 

On the cusp curve 0 = ^ = ^c, the function g{(p) has a point of 
inflection [Fig. 2(b)] and d^g/dtp^, as well as dg/dcp and vanish at 

(p ~ (pp + 2n — arccos[l/(ffp)] = (pc (12c) 

This, in conjimction with t = (v? — ^) /a;, represents the common emission 
time of the three wavefronts that Eire mutually teingential at the cusp curve 
of the envelope. 

In the highly superluminal regime, where f » 1, the sepso-ation of the 
ordinates and 0_ of adjacent maxima and minima in Fig, 2(a) can be 
greater than 27r. A horizontal line <f> = constant will then intersect the curve 
representing g{fp) at more than three points, and so give rise to simtdtane- 
ously received contributions that are made' at 5, 7, • • distinct values of the 
retarded time. In such cases, the sheet 0_ of the envelope (issuing from the 
conical apex of this sxu-face) imdergoes a number of intersections with the 
sheet 0+ before reaching the cusp curve. We shall be concerned in this paper, 
however, mainly with source elements whose distances from the rotation axis 
do not appreciably exceed the radius c/cu of the speed-of-light cylinder and 
so for which the equation g{fp) = (f> has at most three solutions. 

At points of tangency of their fronts, the waves which interfere construc- 
tively to form the envelope propagate normal to the sheets <f) — 4>±{rp^zp) 
of this stirface, in the directions 

h± = {c/u))Vp{<t>± - 4>) 

= erp[fp - fp^(l =F ^h]/R± + e^p/rp + e:,^{zp - z)/R±, (13) 



with the speed c. (e^p , e^^p and e^p are the uni t vectors associated with the 
cylindrical coordinates rp, (pp and zp of the observation point, respectively.) 
Nevertheless, the resulting envelope is a rigidly rotating surface whose shape 
does not change with time: in the (rp, <^p, 2:p)-space, its conical apex is sta- 
tionary at (r, (p, z), and its form and dimensions only depend on the constant 
parameter f . 

The set of waves that superpose coherently to form a particular section 
of the envelope or its cusp, therefore, cannot be the same (i.e. cannot have the 
same emission times) at different observation times. The packet of focused 
waves constituting any given segment of the cusp curve of the envelope, for 
instance, is constantly dispersed and reconstructed out of other waves. This 
one-dimensional caustic would not be unlinaited in its extent, as shown in Fig. 
4, unless the source is infinitely long-lived: only then would the dmration of 
the source encompass the required intervals of emission time for every one of 
its constituent segments. 

II. AMPLITUDES OF THE WAVES GENERATED BY 
A POINT SOURCE 

Our discussion has been restricted so far to the geometrical features 
of the emitted wavefronts. In this section we proceed to find the Lienetrd- 
Wiechert potential for these waves. 

The scalar potential arising from an element of the moving volume 
source we have been considering is given by the retarded solution of the 
wave equation 

V'^Go - d^Go/d(ct')'' = -47rpo, (14a) 

in which 

Po{r\ z\ <') = 5{r' - r)<5(y>' - a;*' - 0)6{z' - z)/r' (146) 

is the density of a point source of unit strength with the trajectory (1). In 
the absence of boundaries, therefore, this potential has the v£due 

C?o(xp, tp) = J d^x'dtVo(x', t')S{tp - <' - |xp - x'|/c)/|xp - x'l (15a) 
dt'S[tp -t'- R{t')/c]/R{t'), (156) 

-oo 



where R{t') is the function defined in (4) (see e.g. Jackson, Classical Electro- 
dynamics, Wiley, New York 1975). 
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If we use (1) to change the integration variable t' in (15b) to v?, and 
express the resulting integrand in terms of the qunatities introduced in (3), 
(5) and (6), we arrive at 

/+00 
dip5[g{ip) - <f>]/R{<p)^ (16) 
•oo 

This can then be rewritten, by formally evaluating the integral, as 

where the angles (pj are the solutions of the transcendental equation g{(p) = 0 
in —oo < (p < +00 and correspond, in conjunction with (1), to the retarded 
times at which the somrce point (r, <p, z) makes its contribution towards the 
value of Go at the observation point (rp^tpp^zp). 

Equation (17) shows, in the light of Fig. 2, that the potential Go of 
a point source is discontinuous on the envelope of the wavefronts: if we 
approach the envelope from outside, the smn in (17) has only a single term 
and yields a finite value for Go, but if we approach this surface from inside, 
two of the 4>jS coalesce at an extremum of g and (17) yields a divergent value 
for Go. Approaching the sheet ^ = 0+ or ^ = <f>^ of the envelope from 
inside this surface corresponds, in Fig. 2, to raising or lowering a horizontal 
line 0 = 00 = const., with 0- < ^o < 0+i until it intersects the cmrve (a) 
of this figure at its maximum or minimimi tangentially. At an observation 
point thus approached, the sum in (17) has three terms, two of which tend 
to infinity. 

On the other hand, approaching a neighboinring observation point just 
outside the sheet <f> — <f>- (say) of the envelope corresponds, in Fig. 2, to 
rsdsing a horizontsJ line 0 = ^o = const., with 4>o < 0-, towards a limiting 
position in which it tends to touch ciurve (a) at its minimum. So long as it 
has not yet reached the limit, such a line intersects curve (a) at one point 
only. The equation g{ip) = 0 therefore has only a single solution = ipoMt iu 
this case which is different from both fp^ and y?- and so at which dgjdtp is 
non-zero (see Fig. 2). The contribution that the source makes when located 
at = (poxxt is received by both observers, but the constructively interfering 
waves that are emitted at the two retarded positions approaching ip^ only 
reach the observer inside the envelope. 

The function Go has an even stronger singularity at the cusp cinrve of 
the envelope. On this cmrve, all three of the (pjS coalesce [Fig. 2(b)] and each 
denominator in the expression in (17) both vanishes and has a vanishing 
derivative {dg/dip = d^g/dip'^ = 0). 
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There is a standard asymptotic technique for evaluating radiation in- 
tegrals with coalescing critical points that describe caustics. By applying 
this technique — which we have outlined in Appendix A — to the integral in 
(16), we can obtain a uniform asymptotic approximation to Go for small 
— <^-|> for points close to the cusp curve of the envelope where Go is 
most singular. The result is 

Go'" 2ci"^(l - X^)~^[pocos(^ arcsinx) - Cigosin(| arcsinx)], \x\ < 1, 

(18) 

and 

Go*'"* - ci'^ix^ - 1)"2 [posinh(|arccosh|x|) 

+ cigosgn(x) sinh(|arccosh|x|)], |xl > li(19) 

where Ci, poi Qo and x are the functions of (r, z) defined in (A2), (A5), (A6) 
and (AlO), and approximated in (A23)-(A30). The superscripts Mn' and 
'out' designate the values of Go inside and outside the envelope, and the 
variable x equals +1 and —1 on the sheets ^ = and <f> = if>- of this 
STurface, respectively. 

The fimction Go°"* is indeterminate but finite on the envelope [c/. 
(A39)], whereas Go''' diverges like y/3ci~^{po =F cigo)/(l - X^)^ as x ±1- 
The singularity structure of Go'" close to the cusp cnrve is explicitly exhibited 
by 

Go- ^ ^(a;/c)(f2f|> - l)-^co^ (z^ - i)^/[co^(fc - 5)' - (^c " <f>y]K (20) 

in which 0 < Zc — z <^ 1, [^c — 01 <^ 1 and 

Co = -^(r^rj, - l)-^(f2,^- l)^(f2 - l)i (21) 

[see (18) and (A22)-(A26)]. It can be seen from this expression that both the 
singularity on the envelope (at which the quantity inside the square brackets 
vanishes) and the singularity at the cusp curve (at which Zc — i and 0c — 0 
vanish) are integrable singularities. 

The potential of a volume source, which is given by the superposition 
of the potentials Go of its constituent volume elements, and so involves in- 
tegrations with respect to (r, ^, z), is therefore finite. Since they are created 
by the coordinated motion of aggregates of p£trticles, the types of sources 
we have been considering cannot, of course, be point-like. It is only in the 
physically imrealizable case where a superluminal source is point-like that its 
potential has the extended singularities described above. 



In fact, not only is the potential of an extended superluminally moving 
source singularity free, but it decays in the far zone like the potential of any 
other source. The following alternative form of the retarded solution to the 
wave equation V^-4o - d^Ao/d{ct)^ = -Airp [which may be obtained from 
(15a) by performing the integration with respect to time]: 

Ao= J d^xp{yi,tp - |x - xp|/c)/|x - xp| (22) 

shows that if the density p of the source is finite and vanishes outside a finite 
volume, then the potential Ao decays Uke |xp|"^ as the distance |xp - x| c:^ 
|xi>| of the observer from the source tends to infinity. 



III. THE BIFURCATION SURFACE OF AN OBSERVER 

Let us now consider an extended somrce which rotates about the z-axis 
with the constant angular frequency a;. The density of such a source — when 
it has a distribution with an unchanging pattern — is given by 

p(r, ^, z, t) = p{r, <p, z), (23) 

where the Lagrangian variable (p is defined by ^ — a;t as in (1), and p can be 
any function of (r, ^, z) that vanishes outside a finite volume. 

If we insert this density in the expression for the reteo-ded scalar poten- 
tial and change the variables of integration from (r, y?, z,t) to (r, ^,z,t), we 
obtain 

^o(xp,tp) = J d^xdtp{y:,t)S{tp - t - |x - xpl/c)/lx - xp| (24a) 
= J rdrd(fidzp{T, <p, z)Go(r, rp, (p - (pp, z - zp), (246) 

where Go is the function defined in (16) which represents the scalar potential 
of a corresponding point source. That the potential of the extended source 
in question is given by the superposition of the potentials of the moving 
source points that consititute it is an advantage that is gained by marking 
the space of source points with the natural coordinates (r, tp, z) of the source 
distribution. This advantage is lost if we use any other coordinates. 

In Sec. II, where the source was point-like, the coordinates (r, ^, z) 
of the source point in Go (r, rp, ^ — ^p, z — zp) were held fixed and we were 
concerned with the behavioxn: of this potential as a function of the coordinates 
(rp^tpp^zp) of the observation point. When we superpose the potentials of 
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the volume elements that constitute an extended source, on the other hand, 
the coordinates {rp^tpp^zp) axe held fixed and we are primarily concerned 
with the behaviour of Go as a function of the integration variables (r, ^, 2). 

Because Go is invariant under the interchange of (r, (p, z) and (rp, ^p, zp) 
if <f> is at the same time changed to —<t> [see (5) and (16)], the singularity of 
Go occTurs on a surface in the (r, ^, z)-space of source points which has the 
same shape as the envelope shown in Fig. 3 but issues from the fixed point 
{rp^(pp^zp) and spirals around the z-aixis in the opposite direction to the 
envelope. [See Fig. 5 in which the light cylinder and the bifurcation surface 
associated with the observation point P are shown for a counterclockwise 
source motion. In this figure, P is located at fp — 9, and only those p8u-ts 
of these surfaces are shown which he within the cyUndrical volume f < 11, 
—1.5 < z — zp < 1.5. The two sheets <f> = <t>±{r^z) of the bifurcation surface 
meet along a cusp (a ciurve of the same shape as that shown in Fig. 4) that 
is tangent to the light cyUnder. For an observation point in the far zone 
(fp 1), the spiralling surface that issues from P undergoes a large number 
of txu-ns — in which its two sheets intersect one another — bofore reaching the 
hght cyUnder.] 

In this paper, we refer to this locus of singularities of Gq as the bifur- 
cation surface of the observation point P. 

Consider an observation point P for which the bifurcation sm^ace in- 
tersects the soinrce distribution, as in Fig. 6. [In Fig. 6, the full curves depict 
the cross section, with the cylinder f = 1.5, of the bifurcation surface of an 
observer located at fp = 3. (The motion of the source is cotmterclockwise.) 
Projection of the cusp cmrve of this bifurcation surface onto the cylinder 
f = 1.5 is shown as a dotted curve, and the region occupied by the soTxrce as 
a dotted area. In this figixre the observer's position is such that one of the 
points (0 = <^C5 ^ = ^Tc) at which the cusp cinrve in question intersects the 
cylinder f = 1.5 — the one with Zc > 0 — is located within the source distri- 
bution. As the radial position rp of the observation point tends to infinity, 
the separation — at a finite distance z^ — z from (0c ? z^ — of the shown cross 
sections decreases like rp~2.] ■ 

The envelope of the wavefronts emanating from a volume element of 
the part of the som"ce that lies within this bifurcation surface encloses the 
point P, but P is exterior to the envelope associated with a source element 
that lies outside the bifurcation surface. 

We have seen that three wavefronts — propagating in diflFerent directions — 
simultaneously pass an observer who is located inside the envelope of the 
waves emanating from a point source, and only one wavefront passes an ob- 
server outside this sTurface. Hence, in contrast to the source elements outside 
the bifmrcation smrface which influence the potential at P at only a single 
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value of the retarded time, this potential receives contributions from each 
of the elements inside the bifurcation surface at three distinct values of the 
retarded time. 

The elements inside but adjacent to the bifurcation surface, for which 
Go diverges, are sources of the constructively interfering waves that not only 
arrive at P simultaneously but also are emitted at the same (retarded) time. 
These somrce elements approach the observer along the radiation direction 
xp — X with the wave speed at the retarded time, i.e. are located at distances 
R{t) from the observer for which 



dR 
dt 



= -c (25) 

t=tp-R/c 



[see (4), (7) and (8)]. Their accelerations at the retarded time 



= T^, (26) 



are positive on the sheet <f> =^ <f>^ of the bifurcation surface and negative on 

The soiu-ce points on the cusp curve of the bifurcation surface, for which 
A = 0 and all three of the contributing retarded times coincide, approach 
the observer — according to (26) — with zero acceleration as well as with the 
wave speed. 

Prom a radiative point of view, the most effective volume elements of 
the superluminal source in question are those that approach the observer 
along the radiation direction with the wave speed and zero acceleration at 
the retarded time, since the ratio of the emission to reception time intervals 
for the waves that are generated by these particular source elements generally 
exceeds unity by several orders of magnitude. On each constituent ring of the 
source distribution that lies outside the light cy Under (r = c/u;) in a plane 
of rotation containing the observation point, there are two volxmae elements 
that approach the observer with the wave speed at the retarded time: one 
whose distance from the observer diminishes with positive acceleration, and 
another for which this acceleration is negative. These two elements are closer 
to one another the smaller the radius of the ring. For the smallest of such 
constituent rings, i.e. for the one that lies on the light cylinder, the two 
volimie elements in question coincide and approach the observer also with 
zero acceleration. 

The other constituent rings of the soiurce distribution (those on the 
planes of rotation which do not pass through the observation point) likewise 
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contain two such elements if their radii are large enough for their velocity 
rcjBy, to have a component along the radiation direction equal to c. On the 
smallest possible ring in each plane, there is again a single volume element — 
at the limiting position of the two coalescing volume elements of the neigh- 
bouring larger rings — that moves towards the observer not only with the 
wave speed but also with zero acceleration. 

For any given observation point P, the efficiently radiating pairs of 
volume elements on various constituent rings of the source distribution col- 
lectively form a surface: the petrt of the bifurcation svirface associated with 
P which intersects the source distribution. The locus of the coincident pairs 
of volume elements, which is tangent to the light cylinder at the point where 
it crosses the plane of rotation containing the observer, constitutes the seg- 
ment of the cusp curve of this bifurcation surface that lies within the source 
distribution. 

Thus the bifurcation surface associated with any given observation 
point divides the volume of the soiu-ce into two sets of elements with dif- 
fering influences on the observed field. As in (18) and (19), the potentials 
Go^ and Gq^^^ of the source elements inside said outside the bifurcation sur- 
face have different forms: the boundary |x(r, rp, (f — ^Pj z — zp)\ = 1 between 
the domains of validity of (18) and (19) delineates the envelope of wavefronts 
when the source point (r,*^, z) is fixed and the coordinates (rp^0p,zp) of 
the observation point are variable, and describes the bifurcation surface when 
the observation point (rp, ^p, zp) is fixed and the coordinates (r, <p, z) of the 
source point sweep a volume. 

The expression (24b) for the scalar potential correspondingly splits into 
the following two terms when the observation point is such that the bifurca- 
tion surface intersects the sorurce distribution: 



where dV = rdrd<pdz, Vin and K>ut designate the portions of the source which 
fall inside and outside the bifurcation surface (see Fig. 6), and Go*" and Go°"* 
denote the different expressions for Go in these two regions. 

Note that the boundaries of the volume VJn depend on the position 
(rp,^p,^p) of the observer: the parameter fp fixes the shape and size of 
the bifurcation surface, and the position (rp, ^p, zp) of the observer specifies 
the location of the conical apex of this surface. When the observation point 
is such that the cusp curve of the bifurcation surface intersects the source 
distribution, the volume Vin is bounded by 0 = <^-, <^ = and the. part of 




(27a) 



(276) 



the source boundary p(r, ^, z) = 0 that falls within the bifurcation surface. 
The corresponding volume Vout is boimded by the same patches of the two 
sheets of the bifiurcation surface and by the remainder of the soxirce boundary. 

In the vicinity of the cusp curve (12), i.e. for - <^ 1 and 0 < Zc - 
z <^1, the cross section of the bifurcation surface with a cylinder f=constant 
is described by 

^^^<t>c c^- (r^ - l)i (4 - \)HfHl - 1)-^ (5c - z) 

± A^f2 _ l)f (^:2^ - l)f (f|,f2 - l)-t(ie - Z)i (28) 

[see (10)-(12) and (A26)]. This cross section, which is shown in Fig. 6, has 
two branches meeting at the intersections of the cusp cm^e with the cylinder 
f — constant whose separation in <^»— at a given (zc — z) — diminishes hke fp^ 
in the limit fp oo. Thus, at finite distances ic — 5 from the cusp curve, 
the two sheets <f> = <f>- and <f> = <f>^ of the bifurcation sxurface coalesce and 
become coincident with the sirrface (f> — §(0- + = C2 as fp — > oo. That 

is to say, the volume Vin vanishes like fp^. 

Because the dominant contributions towards the value of the radiation 
field come from those source elements that approach the observer — along 
the radiation direction — with the wave speed and zero acceleration at the 
retarded time, in what follows, we shall be primarily interested in far-field 
observers the cusp curves of whose bifurcation surfaces intersect the source 
distribution. For such observers, the Green's function hm^ p->oo Go undergoes 
a jump discontinuity across the coalescing sheets of the bifurcation surface: 
the values of x on the sheets ^ = 0±, and hence the functions Go^"*|^=<^_ 
and <Jo**^*U=^+> remain diflferent even in the limit where <f> = <t>- and ^ = 0+ 
coincide [cf. (AlO) and (A39)]. 

IV- DERIVATIVES OF THE RADIATION INTEGRALS AND THEIR 
HADAMARD»S FINITE PARTS 

A. Gradient of the scalar potential 

In this section we begin the calculation of the electric and magnetic 
fields by finding the gradient of the scalar potential Aq, i.e. by calculat- 
ing the derivatives of the integral in (27a) with respect to the coordinates 
(rp,¥?p,zp) of the observation point. 

If we regard its singular kernel C?o as a classical function, then the inte- 
gral in (27a) is improper and caimot be diflFerentiated under the integral sign 
without characterizing and duly handling the singularities of its integrand. 
On the other hand, if we regard Go as a generalized function, then it would 
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be mathematically permissible to interchange the orders of differentiation 
and integration when calculating VpAo, 

This interchange results in a new kernel VpCo whose singularities are 
non-integrable. However, the theory of generalized functions prescribes a 
well-defined procedure for obtaining the physically relevant value of the re- 
sulting divergent integral, a procedure involving integration by parts which 
extracts the so-called Hadamard^s finite part of this integral (see e.g. Hoskins, 
Generalised Functions, Ellis Horwood, London 1979). Hadamaxd's finite part 
of the divergent integral representing VpAo yields the vahxe that we would 
have obtained if we had first evaluated the original integreJ for Aq as an 
explicit function of {rp,(pp^ zp) and then differentiated it. 

From the standpoint of the theory of generalized functions, therefore, 
differentiation of (27a) yields 

Vp^o = J dVpVpGo = (Vp^o)in + (Vp^o)out, (29a) 

in which 

(Vp^o)in.out = / dypVpGo'"'°"*. (296) 

Since p vanishes outside a finite volume, the integral in (27a) extends over all 
values of (r, ^, z) Eind so there is no contribution from the limits of integration 
towstrds the derivative of this integral. 

The kernels VpGo*"'^"* of the above integrals may be obtained from 
(16). Applying Vp to the right-hand side of (16) and interchanging the 
orders of differentiation and integration, we obtain an integral representation 
of VpGo consisting of two terms: one arising fi:om the differentiation of R 
which decays like rp~^asrp— >oo and so makes no contribution to the field 
in the radiation zone, and another that arises from the differentiation of the 
Dirac delta function and decays less rapidly than rp"^. For an observation 
point in the radiation zone, we may discard terms of the order of rp"^ and 
write 

/H-oo 
d<pR-^S'{g - 0)n, fp » 1, (30) 
■oo 

in which 5' is the derivative of the Dirac delta fimction with respect to its 
argximent and 

n = erp [fp - f cos{(p - (pp)]/R + e^^/fp + e^p (zp - z)/R. (31) 

Equation (30) yields V/>Go*" or VpGo^"* depending on whether 0 lies within 
the interval (^_, 0+) or outside it. 
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If we now insert (30) in (29b) and perform the integrations with respect 
to by parts, we find that 

(Vp>lo)in - Wc) J d(f>dp/dv>Gi'''], fp » 1, 

* (32) 

£Lnd 

(Vp^o)out ^ (uj/c) rdTdz[ [pGi^"*] 

+ {f ' + P^)d<f>dp/d<pGx''''^]. fp » 1, (33) 

in which S stands for the projection of Vin onto the (r,z)-plane, and Gi'" 
and Gi°"*^ are given by the values of 

/ + 00 
dcpR-^5{g -<f>)n= ^'^l^d/M'^^ (34) 
-oo ^^^^ 

for <f> inside and outside the interval (^- , i^^), respectively. 

Like Go''', the Green's function Gi''' diverges on the bifurcation surface 
= ^i, where dg/d(p vanishes, but this singularity of Go*'* is integrable 
so that the value of the second integral in (32) is finite (see Sec. II and 
Appendix A). Hadamard's finite part of (Vp>lo)in (denoted by the prefix 
Fp) is obtaind by simply discarding those 'integrated' or boundary terms in 
(32) which diverge. Hence, the physically relevant quantity Fp{(VpAo)in} 
consists— in the far zone — of the volume integral in (32). 

Let us choose an observation point for which the cusp curve of the bi- 
furcation smrface intersects the source distribution (see Fig. 6). When the 
dimensions L) of the source are negligibly smaller than those of the bifur- 
cation surface (i.e. when L <^rp and so Zc — z <^ rp throughout the source 
distribution) the functions 01'"'°"** in (32) and (33) can be approximated by 
their asymptotic values (A34) and (A35) in the vicinity of the cusp curve 
(see Appendix A). 

According to (A34), (A36) and (A44), Gi''' decays like pi/ci^ = 0(1) 
at points interior to the bifurcation surface where limiz^^-^cxs X remains finite. 
Since the separation of the two sheets of the bifurcation surface diminishes 

like fp^ within the source [see (28)], it therefore follows that the volxmae 

_ s 

integral in (32) is of the order of 1 x fp^ , a result which can also be inferred 
from the far-field version of (A34) by explicit integration. Hence, 

Fp{(Vp>lo)in} = 0(f;^), fp » 1, (35) 



16 



decays too rapidly to make any contribution towards the value of the electric 
field in the radiation zone. 

Because Gi^"^ is, in contrast to Gi^'*, finite on the bifurcation surface, 
both the surface and the volume integrals on the right-hand side of (33) have 
finite values. Each component of the second term has the same structiu^e as 
the expression for the potential itself and so decays like rp~^ (see the ultimate 
paragraph of Sec. II). But the first term — which would have cancelled the 
correspoding boundary term in (32) and so would not have survived in the 
expression for Vp^o had the Green's function Gi been continuous — behaves 
differently firom any conventional contribution to a radiation field. 

Insertion of (A39) in (33) yields the following expression for the asymp- 
totic value of this boundary term in the limit where the observer is located 
in the fgtr zone and the source is localized about the cusp curve of his (her) 
bifurcation sm~face: 



rdr(iz[pGi°"*]^"' ~ |ci-2/rdrdz[pi(pU^ -/£)U_) + 2ciqi(p|^++H^_)]. 

(36) 

In this limit, the two sheets of the bifurcation sxrrface are essentially coinci- 
dent throughout the domain of integration in (36) [see (28)]. So the difference 
between the values of the source density on these two sheets of the bifurcation 
surface is negligibly small (^ fp^) for a smoothly distributed source and the 
functions p\^^ appearing in the integrand of (36) may correspondingly be 
approximated by their conunon limiting value Pbs('">^) on these coalescing 
sheets. 

Once the functions pj^^. are approximated by Pbs(^> ^) aiid qi by (A41), 
equation (36) yields an expression which can be written, to within the leading 
order in the far-field approximation fp » 1 [see (A44) and (A45)], as 

J rdrdz[pGi^"^]^'' - 2§(c/a;)2fp^ J^^ df{f^ - l)-^iii 

J Zc — Xff a;/c 

-2t(c/a;)2f;» r^df(f2 -l)-ini(Wc)»(/>b8), (37) 

with ^ 

(pbs)(r)= / dr}pbs{T,z)\^^^^_^,j^^, (38) 

J 0 

where Zc" Lz{t) < z < z^ and r< < r < r> are the intervals over which 
the bifurcation surface intersects the source distribution (see Fig. 6). The 
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quantity (/c?bs)('") ™ay be interpreted, at any given r, as a weighted average — 
over the intersection of the coalescing sheets of the bifxu-cation surface with 
the plane z = Zc — V^L^ — of the source density p. 

The right-hand side of (37) decays like rp~t as rp oo. The sec- 
ond term in (33) thus dominates the first term in this equation, and so the 
quantity (Vpj4o)out itself decays like Tp~^ in the far zone. 

B. Time derivative of the vector potential 

Inasmuch as the charge density (23) has an unchanging distribution 
pattern in the (r, ^, 2:)-frame, the electric current density associated with the 
moving source we have been considering is given by 

j(x, t) = rLjp{r, <p, z)ey, (39) 

in which rcoBtp = ru;[— sin(¥? — Vp)^rp + cos(v? — Vp)^ipp\ is the velocity of 
the element of the som"ce pattern that is located at {r^(p^z). This current 
satisfies the contintiity equation dp/d{ct) -t- V • j = 0 automatically. 

In the Lorentz gauge, the retarded vector potential corresponding to 
(24a) has the form 

A(xp,<p) = c-^ j d^xdtj(yi, t)6{tp - < - |x - xp|/c)/|x - xp|. (40) 

If we insert (39) in (40) and change the variables of integration firom (r, ^, z, t) 
to (r, y?, z,(p)j as in (24), we obtain 

A = y^ dVfp{r,(p,z)G2{r,rp,i}> - (pp,z~ zp), (41) 

in whcih dV = rdrdcpdz^ the vector G2 — which plays the role of a Green's 
function — is given by 

dv>e^5[g(ip) - <f>yR{<p) = Yl R-^\dg/d>fi\-^%^, (42) 

-OO 

and g and (pjS are the same quantities as those appearing in (17) (see also 
Fig. 2). 

Because (17), (34) and (42) have the factor \dg/d(p\~^ in common, the 
fxmction G2 has the same singularity structure as those of Go and Gi: it 
diverges on the bifurcation surface dg/d<p = 0 if this surface is approached 
from inside, and it is most singular on the cusp cuirve of the bifurcation surface 
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where in addition d^g/dip^ = 0. It is, moreover, described by two different 
expressions, and 62°^*, inside and outside the bifurcation srurface whose 
asymptotic values in the neighbourhood of the cusp curve have exactly the 
same functional forms as those found in (18) and (19); the only difference 
being that po and go in these expressions are replaced by the P2 and q2 given 
in (A37) (see Appendix A). 

As in (29), therefore, the time derivative of the vector potential has the 
form d A/dtp = {dA/dtp)in + (dA/dtp)out with 

(aA/a<p)in,out = -u; I dVfpaGs'"'^"' (43) 

when the observation point is such that the bifurcation surface intersects the 
source distribution. 

The functions G2^"'°"* depend on <pp and (p in the combination (p — 0p 
only. We can therefore replace djdCpp in (43) by —d/d(p and perform the 
integration with respect to (p by p£ui;s to arrive at 

{dA/dtp)in ^cj^ drdzf^ { [pG2^"] - J^'^ d<f>dp/d<pG2''' } , (44) 

and 

{dA/dtp)ont ^-^/^ { [PG2^"*] 

+ J^^)d<f>dp/dv>G2''''^]^ (45) 

For the same reasons as those given in the psiragraphs following (32) and 
(33), Hadamard's finite part of {dA/dtp)in consists of the volume integral 
in (44) and is of the order of fp^ [note that according to (A37) and (A42), 
P2 > CiQa and P2/ci^ = 0(1)]. The volume integral in (45), moreover, 
decays like fp"^, as does its coimterpart in (33). 

The part of dA/dtp that decays more slowly than conventional contri- 
butions to a radiation field is the boundary term in (45). The asymptotic 
value of this term is given by an expression similar to that appearing in (36), 
except that pi and qi are replaced by p2 and q2. Once the quantities p|^^ 
and q2 in the expression in question are approximated by pbs and by (A42), 
as before, it follows that . 

{dA/dtp)out c J^drdzf^^pG2''''^Yj |c drdzf2pbsCi"^q2 

- -4i<^yu;)fpe^, i?; drrHf^ - 1)"^ S^^L.^f^dzf^o " i)-^Pbs<46) 
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This behaves like fp^ as fp cxd since the i-quadrature in (46) has the 
finite value 2{LzU)/c)'2 (p^^) in this limit [see (37) et seg.]. 
Hence, the electric field vector of the radiation 

E = -Vp^o - dA/d{ctp) - -c-\dA/dtp)out 

~ 4(c/u;)f;'e^, jT> dff2(f2 - irHL,u;/c)Hp^s) (47) 

itself decays like rp"5 in the far zone: as we have aheady seen in Sec. IV (A), 
the term Vp-Ao heis the conventional rate of decay rp~^ and so is negligible 
relative to (^A/5ip)out. 



C. Curl of the vector potential 

There are no contributions from the limits of integration towards the 
curl of the integral in (41) because p vanishes outside a finite volimie and 
so the integrsJ in this equation extends over all values of (r, ^, z). Hence, 
differentiation of (41) yields 

B = Vpx A = Bin + Bout, (48a) 

in which 

Bin,out = / dVfpVp X G2*"*°"*. (486) 

Operating with VpX on the first member of (42) and ignoring the term 
that decays like rp"^, as in (30), we find that the kernels Vp x 62^^^ and 
Vp X G2*^"* of (48b) are given — in the radiation zone — ^by the values of 

/+00 
dipR'^5'{g - (/>)& X e^, fp » 1, (49) 
■00 

for (f> inside and outside the interval (0_,0-|.), respectively, [n is the xmit 
vector defined in (31).] 

Insertion of (49) in (48) now yields expressions whose ^quadratures 
can be evaluated by parts to arrive at 

Bin 



and 



f drdzf^{ ~ [pG3H^_^'" + f ^ d<t>dp/d(pG3^}, rp » 1, (50) 

+ (y^^'+^^')d09p/a^G3*'"*}, fp»l, (51) 
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where and Ga""^ stand for the values of 

-00 ^^^^ 

inside and outside the bifurcation surface. 

Once again, owing to the presence of the factor \dg/d<p\^'^ in Gs**^, 
the first term in (50) is divergent so that the Hadamard's finite part of Bin 
consists of the volume integral in this equation, an integral whose magnitude 
is of the order of fp^ [see the paragraph containing (35) and note that, 
accroding to (ASS) and (A44), pa » ciQa and pa/ci^ = 0(1)]. The second 
term in (51) has— like those in (33) and (45)— the conventional rate of decay 
fp^. Moreover, the surface integral in (51) — which would have had the same 
magnitude as the surface integral in (50) and so would have cancelled out 
of the expression for B had G3'" and Gs"""* matched smoothly across the 
bifurcation surface — decays as slowly as the corresponding term in (45). 

The asymptotic value of G3 for source points close to the cusp curve of 
the bifurcation sinrface has been calculated in Appendix A. It follows from 
this value of G3 and fi^om (51), (52), (A40), (A44) and (A45) that, in the 
radiation zone, 

J drdzf"" [pG3^"*] ^ I /s drdzf2pbsCi-iq3 
^ ^-i^M^p' r> dfr^f^ - 1)-^ f!:^L.u,/cd^i^c - i)-^Pbsn3 (53) 

to within the order of the approximation entering (37) and (46)- 

The far-field version of the radial unit vector defined in (31) assinmes 

the form , . 

Hm n| . = f-^Br^ - (1 - f-2)le., (54) 

on the cusp cturve of the bifmrcation sinrface [see (12b), (13) and (A27), and 
note that the position of the observer is here assimied to be such that the 
segment of the cusp curve lying within the source distribution is described 
by the expression with the plus sign in (12b), as in Fig. 6]. So, ns equals 
fi X e^p in the regime of validity of (53) [see (A45)]. Moreover, n can be 
replaced by its far-field value 

hc^{rpBrj>+zpezp)/Rpy Rp oo, (55) 

if it is borne in mind that (53) holds true only for an observer the cusp curve 
of whose bifurcation surface intersects the source distribution. 
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Once xiz in (53) is approximated by n x Bi^^ and the resiilting z- 
quadrature is expressed in terms of (pbs) [see (38)], this equation reduces 
to 

B n X E, (56) 

where E is the electric field vector earlier found in (47) . Equations (47) and 
(56) jointly describe a radiation field whose polarization vector lies along the 
direction of motion of the source, 6^,^ . 

Note that there has been no contribution toward the values of E and 
B from inside the bifxurcation surface. These quantities have arisen in the 
above cgJculation solely firom the jump discontinuities in the values of the 
Green's functions Gi*""*, and across the coalescing sheets of 

the bifurcation surface. We would have obtained the same results had we 
simply excised the vanishingly small volume limrp-^00 from the domains 
of integration in (29), (43) and (48). 

Note also that the way in which the familiar relation (56) has emerged 
from the present analysis is altogether different from that in which it ap- 
peeo^ in conventional radiation theory. Essential though it is to the physicsil 
requirement that the directions of propagation of the waves and of their en- 
ergy should be the same, (56) expresses a relationship between fields that are 
here given by non-spherically decaying svirface integrals rather than by the 
conventional volume integrals that decay like rp~^. 

V. A PHYSICAL DESCRIPTION OF THE EMISSION PROCESS 

Expressions (47) and (56) for the electric and magnetic fields of the 
radiation that arises from a charge-current density with the components (23) 
and (39) imply the following Poynting vector: 

S ~ i^7r-ic(c/a;)2fp^[ jT> dff^r^ - {L,uj/c)Hp^,)]^n. (57) 

In contrast, the magnitude of the Poynting vector for the coherent cyclotron 
radiation that would be generated by a macroscopic lump of charge, if it 
moved subluminally with a centripetal acceleration cu;, is of the order of 
{{p}L^)^u^/(cRj>^) according to the L£u:mor formula, where represents 
the volume of the source and (p) its average charge density. The intensity 
of the present emission is therefore greater than that of even a coherent 
conventional radiation by a factor of the order of {L2/L){Li*)/c)~^{Rp/L)^ a 
factor that ranges from 10^^ to 10^^ in the case of pulsars for instance. 

The reason this ratio has so large a value in the far field (Rp/L ^ 1) 
is that the radiative characteristics of a volume-distributed som"ce which 
moves faster than the waves it emits are radically different from those of a 
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corresponding source that moves more slowly than the waves it emits. There 
axe source elements in the former case that approach the observer along the 
radiation direction with the wave speed at the retarded time. These lie on 
the intersection of the source distribution with what we have here called the 
bifurca-tion surface of the observer (see Figs. 5 and 6): a surface issuing from 
the position of the observer which has the same shape as the envelope of 
the wavefronts emanating from a soxarce element (Figs. 1 and 3) but which 
spirals around the rotation axis in the opposite direction to this envelope 
and resides in the space of source points instead of the space of observation 
points. 

The source elements inside the bifiurcation smrface of an observer make 
their contributions towards the observed field at three distinct instants of 
the retarded time. The values of two of these retarded times coincide for an 
interior source element that lies next to the bifurcation surface. This limit- 
ing value of the coincident retarded times represents the instant at which the 
component of the velocity of the source point in question equals the wave 
speed c in the direction of the observer. The third retarded time at which a 
soiu*ce point adjacent to — just inside — the bifurcation surface makes a con- 
tribution is the same as the single retarded time at which its neighbouring 
source element just outside the bifurcation surface makes its contribution 
towards the observed field. (The soiurce elements outside the bifurcation sxir- 
face make their contributions at only a single instant of the retarded time). 

At the instant marked by this third value of the retarded time, the 
two neighbouring source elements — just interior and just exterior to the bi- 
furcation surface — have the same velocity, but a velocity whose component 
along the radiation direction is different fi:om c. The velocities of these two 
neighbouring elements are, of course, equal at any time. However, at the 
time they approach the observer with the wave speed, the element inside 
the bifurcation surface makes a contribution towards the observed field while 
the one outside this surface does not: the observer is located just inside the 
envelope of the wavefronts that emanate from the interior source element but 
just outside the envelope of the wavefronts that emanate from the exterior 
one. ThTis, the constructive interference of the waves that are emitted by 
the source element just outside the bifurcation smrface takes place along a 
caustic which at no point propagates past the observer at the conical apex 
of the bifurcation sxurface in question. 

On the other hand, the radiation effectiveness of a source element which 
approaches the observer with the wave speed at the retarded time is much 
greater than that of a neighboxaring element the component of whose velocity 
along the radiation direction is subluminal or superluminal at this time. This 
is because the piling up of the emitted wavefronts along the line joining the 
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source and the observer makes the ratio of emission to reception time intervals 
for the contributions of the lumineJly moving source elements by many orders 
of magnitude greater than that for the contributions of any other elements. 
As a result, the radiation effectiveness of the various constituent elements 
of the source (i.e. the Green's function for the emission process) undergoes 
a discontinuity across the boundeiry set by the bifiu-cation surface of the 
observer. 

The integral representing the superposition of the contributions of the 
various volume elements of the source to the potential thus entails a discon- 
tinuous integrand. When this volume integral is differentiated to obtain the 
field, the discontinuity in question gives rise to a boundary contribution in 
the form of a surface integral over its locus. This integral receives contri- 
butions from opposite faces of each sheet of the bifurcation surface which 
do not cancel one smother. Moreover, the contributions arising from the 
exterior faces of the two sheets of the bifurcation siurface do not have the 
same value even in the limit i2p — > oo where this sin-face is infinitely large 
and so its two sheets are — throughout a localized soiirce that intersects the 
cusp — coalescent. Thus the resulting expression for the field in the radiation 
zone entails a sm-face integral such as that which would arise if the source 
were two-dimensional, i.e. if the source were concentrated into an infinitely 
thin sheet that coincided with the intersection of the coalescing sheets of the 
bifmrcation surface with the source distribution. 

For a two-dimensional source of this type — whether it be real or a vir- 
tual one whose field is described by a surface integral— the near zone (the 
Presnel regime) of the radiation can extend to infinity, so that the amplitudes 
of the emitted waves are not necesssirily subject to the spherical spreading 
that normally occurs in the far zone (the Praunhofer regime). The Presnel 
distande which marks the boundary between these two zones is given by 
Rp ^ L±^/L\\, in which L± and are the dimensions of the source perpen- 
dicular and parallel to the radiation direction. If the source is distributed 
over a sin^face and so has a dimension L|| that is vanishingly small, therefore, 
the Presnel distance Rr tends to infinity. 

In the present case, the surface integred which arises from the disconti- 
nuity in the radiation effectiveness of the source elements across the bifurca- 
tion siurface has an integrand that is in txim singular on the cusp curve of this 
surface. This has to do with the fact that the source elements on the cusp 
curve of the bifurcation surface approach the observer along the radiation 
direction not only with the wave speed but also with zero acceleration. The 
ratio of the emission to reception time intervals for the signals generated by 
these elements is by several orders of magnitude greater even than that for 
the elements on the bifin-cation surface. When the contributions of these 
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elements sure included in the surface integral in question, i.e. when the obser- 
vation point is such that the cusp curve of the bifurcation surface intersects 
the source distrbution (as shown in Fig. 6), the value of the resulting im- 
proper integral turns out to have the dependence i2p"^, rather than i2p"^, 
on the distance Rp of the observer from the source. 

This non-spherically decaying component of the radiation is in addi- 
tion to the conventional component that is conciu*rently generated by the 
remaining volume elements of the soiurce. It is detectable only at those ob- 
servation points the cusp ciurves of whose bifvircation smrfaces intersect the 
source distribution. It appears, therefore, as a spiral-shaped wave packet 
with the same azimuthal width as the ^extent of the sorirce. For a som*ce 
distribution whose superluminal portion extends from f = ltof = f>>l, 
this wave packet is detectable — by an observer at infinity — within the angles 
^TT — arccosf>^ < 0p < + arccosf>^ from the rotation axis: projection 
(12b) of the cusp cmrve of the bifurcation siurface onto the (r, z)-plane reduces 
to cot^p = (f^ - 1)= in the limit Rp 00, where Op = arctan(rp/zp) [also 
see (54)]. 

Because it comprises a collection of the spiralling cusps of the envelopes 
of the wavefronts that are emitted by varioxis source elements, this wave 
packet has a cross section with the plane of rotation whose extent and shape 
match those of the source distribution. It is a diffraction-free propagating 
caustic that — when detected by a far-field observer — would appear as a pulse 
of duration A^/a;, where Atp is the azimuthal extent of the soiurce. 

Note that the waves that interfere constructively to form each cusp, 
and hence the observed pulse, are different at different observation times: 
the constituent waves propagate in the radiation direction n with the speed 
c, whereas the propagating caustic that is observed, i.e. the segment of the 
cvisp curve that p£isses through the observation point at the observation time, 
propagates in the azimuthal direction e^^ with the phase speed rpuj. 

The fact that the intensity of the pulse decays more slowly than pre- 
dicted by the inverse square law is not therefore incompatible with the con- 
servation of energy, for it is not the same wave packet that is observed at 
different distances from the source: the wave packet in question is constantly 
dispersed and re-costructedted out of other waves. The cusp ctirve of the 
envelope of the wavefronts emanating from an infinitely long-lived source is 
detectable in the radiation zone not because any segment of this curve can be 
identified with a caustic that has formed at the source and has subsequently 
travelled as an isolated wavepacket to the radiation zone, but because certain 
set of waves superpose coherently only at infinity. 

Relative phases of the set of waves that are emitted dinring a limited 
time interval is such that these waves do not, in general, interfere construc- 
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tively to form a cusped envelope until they have propagated some distance 
away from the source. The period in which this set of waves has a cxisped 
envelope and so is detectable as a periodic train of non-spherically decaying 
pulses, would of course have a limited duration if the source is short-lived. 

Thus, pulses of focused waves may be generated by the present emis- 
sion process which not only are stronger in the far field than any previously 
studied class of signals, but which can in addition be beamed at only a select 
set of observers for a limited interval of time. 

VI. DESCRIPTION OF EXAMPLES OF THE APPARATUS 

An apparatus can be designed for generating such pulses, in accordance 
with the above theory, which basically entails the simple components shown 
in Figs. 7(a) and 7(b). 

Referring to the example of Fig. 7(a), a linear dielectric rod 1 of length 
I is provided with an ctrray of electrodes 2, 3 arranged opposite one another 
along its length with n/t electrodes per unit length. In use, a voltage po- 
tential is applied across the dielectric rod 1 by the electrodes 2, 3, with each 
pair of electrodes 2, 3, in the array being activated in tvim to generate a 
polarisation region with the fronts 5. By rapid application and removal of a 
potential voltage to electrodes 2, 3, this polarised region can be set in motion 
with a superlxuninal velocity. Creating a voltage across a pair of electrodes 
polarises the material in the rod between the electrodes. The electrodes can 
be controlled independently, so that the distribution pattern of polarisation 
of the rod as a function of length edong the rod is controlled. 

By varying the voltage across the electrode pairs as a function of time, 
this polarisation pattern is set in motion. For example, neighbouring elec- 
trode pairs can be turned on with a time interval of At between them, start- 
ing from one end of the rod. Thus, at a snapshot in time, part of the rod 
is polarised (that part lying between electrode pairs with a voltage across 
them) and part of it is not polarised (that part lying between electrode pairs 
without a voltage across them). These regions are separated by *'polarisa- 
tion fronts" which move with a speed of l/(nAt). With suitable choices of n 
and A* the polarisation fronts can be made to move at any speed (including 
speeds faster than the speed of light in vacuo). The polarisation fronts can 
be accelerated through the speed of light by changing At with time- 
Figure 7(b) shows another example of the invention, the one ansJysed 
above. In this example, the dielectric rod is formed in the shape of a ring. 
Figxrre 7(b) is a plan view showing electrodes 2, and has electrodes 3 disposed 
below the rod 1. For a ring of radiiis r and a polarisation pattern that moves 
aroxmd the ring with an angiilar frequency a;, the velocity of the charged 
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region is ru). In this example, ruj is greater than the speed of light c so that 
the moving polarisation pattern emits the radiation described with reference 
to Figures 1 to 6. 

The voltages across neighboin-ing electrode pairs have the same time 
dependence (their period is 27r/a;) but, as in the rectilinear case, there is 
a time difference of Ai between them. The polarisation pattern must move 
coherently siroimd the ring, i.e. must move rigidly with an unchanging shape; 
this would be the case if n At = 2'kN/uj, where n is the number of electrodes 
around the ring and N an integer. Within the confines of this condition, the 
time dependence of the voltage across each pair of electrodes can be chosen 
at will. The exact form of the adopted time dependence would allow, for 
example, the generation of harmonic content and structure in the sotirce. 

The electrodes are driven by an atrray of similar oscillators, an array in 
which the phase difference between successive oscillators has a fixed value. 
There are several ways of implementing this: 

a single oscillator may be used to drive each electrode through progres- 
sively longer delay lines; 

each electrode pair may be driven by an individual oscillator in an array 
of phgLse-locked oscillators; or 

the electrode pairs may be connected to points aroxmd a circle of radius 
r which lies within — and is coplanar with — an annular waveguide, a waveg- 
uide whose normal modes include an electromagnetic wave train that prop- 
agates longitudinally around the circle with an angular frequency lj > c/r. 

For a dielectric rod in the shape of a ring of diameter 1 m, oscillators 
operating at a fi*equency of 100 MHz would generate a superltmainally moving 
polarisation pattern. The required oscillator frequencies oxe esisily obtainable 
using standard laboratory equipment, and any material with an appreciable 
polarizability at MHz frequencies would do for the medium. 

With oscillators operating at frequeiicies of 1 GHz (also available), the 
size of the device would be about 10 cm across; applications demanding 
portability are therefore viable. 

Unhke in conventional emissions, the spectrum of the Green's function 
for the present radiation process has no exponential cut-off at high frequen- 
cies. The high-frequency cut-off of the spectrxun of this new t3T>e of radiation 
is determined solely by that of the spectrum of the electric current density j. 
Our ability to control the scale of variations of j in ^ provides, therefore, the 
possibility of using the apparatus described above to generate emission over a 
chosen frequency range from the radio to ultraviolet or even x-ray. (This ap- 
phes also to the radiation arising from the subluminal part 0 < 1 — ra;/c <§C 1 
of the rotating source when u) is such that the speed-of-light cylinder r = c/w 
falls within the dielectric ring of Fig. 7b.) 
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APPENOIX A: ASVMPTOTIC EXPANSIONS OF 
THE GREEN'S FUNCTIONS 

In this Appendix, we calculate the leading terms in the asymptotic 
expansions of the integrals (16), (34), (42) and (52) for small ^+ - i.e. 
for points close to the cusp curve (12) of the bifurcation surface (or of the 
envelope of the wavefronts). The method — originally due to Chester et al 
(Proc. Camb. Phil. Soc, 54, 599, 1957) — which we use is a standard one 
that has been specifically developed for the evaluation of radiation integrals 
involving caustics (see Ludwig, Comm. Pure Appl. Maths, 19, 215, 1966). 
The integrals evaluated below all have a phase function g{ip) whose extrema 
((p = ipj^) coalesce at the caustic (12). 

As long as the observation point does not coincide with the soxirce 
point, the function g{ip) is analytic and the following transformation of the 
integration variables in (16) is permissible: 

9{<p) = - ci'^i^ + C2, {Al) 

where u is the new variable of integration and the coeflBcients 

ci = (|)i(<^+-i?i.)^ and C2 = ^((^+ + 0-) {A2) 

are chosen such that the values of the two functions on opposite sides of (Al) 
coincide at their extrema. Thus an alternative exact expression for Go is 

/H-oo 
dufo{y)5{\u^ - cx^v + C2 - 0), (^3) 
-OO 

in which 

fo{y)^R'^dip/dv. (A4) 

Close to the cusp curve (12), at which ci vanishes and the extrema u = 
±ci of the above cubic fxmction are coincident, fo{^) may be approximated 
by po H- Qo^^ with 

Po=|(/oU,+/oL=_,,), (A5) 

and 

go=icx-i(/oL=,^.-/oL=_,,). (^6) 
The resulting expression 

/+00 
dv{po + qo^)5{\v^ - ci^i/ + C2 - 0) (A7) 
-co 
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will then constitute, according to the general theory, the leading term in the 
asymptotic expansion of Go for small Ci. 

To evaluate the integral in (A7), we need to know the roots of the cubic 
equation that follows from the vanishing of the argument of the Dirac delta 
function in this expression- Depending on whether the observation point is 
located inside or outside the bifiurcation surface (the envelope), the roots of 

-Ci^Z/ + C2 = 0 (^8) 

are given by 

V = 2ci cos(|n7r + | arccosx), |x| < 1? (-^9^) 
for n =0, 1 and 2, or by 

u = 2cisgn(x) cosh(|arccosh|x:|), Ix! > 1, (^9^) 
respectively, where 

X = - + <t>-)]/[W^ - <!>-)] = - c2)/ci^ (Aio) 

Note that x equals +1 on the sheet <^ = <^^_ of the bifm-cation surface (the 
envelope) and —1 on <f> = 

The integral in (A7), therefore, has the following value when the obser- 
vation point lies inside the bifurcation surface (the envelope): 



dv5{\u^-cr''u+C2) = ELoCi-^'Ncos^dnTT+i arccosx)-l| , Ixl < 1- 

(^11) 

Using the trignometric identity 4 cos^ a — 1 = sin 3a/ sin a, we can write this 
as 

+00 ^ 2^ 

du5{\v^ - Ci^i^H- C2) = ci"^(l - x^)"= zZ I sin(|n7r + | arccosx)| 

ns=0 



= 2ci-2(l - x")-"^ cos(i arcsinx), \x\ < 1, (-A12) 



in which we have evaluated the sum by adding the sine functions two at a 
time. 

When the observation point lies outside the bifurcation surface (the 
envelope), the above integral receives a contribution only from the single 
value of u given in (A9b) and we obtain 



+00 

du5{\u^ - oi'^v + C2) = Ci-2(;^2 _ jj-i sinh(|arccosh|x|), \x\ > 1, 

(^13) 



/+OC 
-00 



29 



where this time we have used the identity 4cosh^a — 1 = sinhSa/ sinha. 

The second part of the integral in {A7) can be evaluated in exactly the 
same way. It has the value 



duuSilu^ " ci^u + C2) = 2ci-^(l - x^)~^ I sm(§n7r + | arccosx)| 

n=0 

X cos(|n7r + 5 arccosx) 
= -2ci-ni - X^r^ sin(| arcsinx), Ixl < 1, (^14) 



when the observation point lies inside the bifurcation surface (the envelope), 
ajid the value 



4-00 

diyuSi^u^-ci^u+cz) = ci-Hx:^-l)"*sgn(x)sinh(|arccosh|x|),|xl > 1, 

(^15) 

when the observation point lies outside the bifurcation surface (the envelope). 

Inserting (A12)-(A15) in (A7), and denoting the values of Go inside 
and outside the bifurcation surface (the envelope) by Go""* and Go°"*, we 
obtain 

Go'" 2ci"^(l - X^)"^bocos(| arcsinx) - cigosin(f arcsinx)], Ixl < 1, 

(^16) 

and 

Go"""* ci'^ix^ - l)-2 [posinh(iarccoshlxl) + cigosgn(x) sinh(|arccosh|xl)], 

|X|>1, (^17) 

for the leading terms in the asymptotic approximation to Gq for small ci. 

The function fo{^) iii terms of which the coefficients po and go are 
defined is indeterminate at i/ = ci and u = -cit differentiation of (Al) yields 
d<p/du = (i/^ — ci^)/{dg/dip) the zeros of whose denominator at = and 
<p = (p^ respectively coincide with those of its numerator at u — +ci and 
J/ = — ci. This indeterminacy can be removed by means of THopital's rule 
by noting that 



dtp 
1^ 



I/=±Cl 



dg/dip 



2v 



(^18) 



i.e. that 

dip 
dv 
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in which we have calculated {d^g/dip^)^^^ from (7) and (8). The right-hand 
side of (A19) is, in turn, indeterminate on the cusp curve of the bifurcation 
surface (the envelope) where ci = A = 0. Removing this indeterminacy by 
expanding the numerator in this expression in powers of A4, we find that 
d(p/diy assumes the value 23 at the cusp c\n:ve. 

Hence, the coefficients po and go tiiat appear in the expressions (A8) 
and (A9) for Go are explicitly given by 

PO = (u;/c)(ici)i(J?:^ +^;^)A-i, {A20)) 

and J 

90 = {uj/c){2c^rHRZ' - ^;^)A-i {A21) 

[see (A4)-(A6) and (A19)]. 

In the regime of validity of (AS) and (A9), where A is much smaller 
thein (fpf^ ~ 1)^> leading terms in the expressions for R±y ci^ po and go 
are 

R± = (f|>f2 - 1)^ ± (f|>f2 - 1)-^ A^ + 0(A), {A22) 

ci = 2-3(f|>f2 - I)-^aWO(A), {A23) 

PO. = 23 (a;/c)(f 2f|, - 1)-^ + O(A^), (^24) 

and 

go = 2-i(a;/c)(f2f|> - 1)"^ + O(A^). {A25) 

These may be obtained by using (9) to express z everywhere in (10), (11) 
and (A2) in terms of A and f, and expanding the resulting expressions in 
powers of A2. The quantity A in turn has the following value at points 
0 < ic - i «: (^P - - 1) = : 

A = 2{f% - l)Hf2 - l)i{z^ - z) + 0[(zc - z)% {A26) 

in which ic is given by the expression with the plus sign in (12b). 

For an observation point in the far zone (fp » 1), the above expressions 
reduce to 

R^ c^ ffp, ci - 26 (ffj>)- = (1 - f-2)^(5c - 5)^, {A27) 

A - 2fp(f2 - l)i (ie - 5), (^28) 
Po ^ 23(a;/c)(fpf)-^ go ^ 2^3 (a;/c)(fpf)-^ (^29) 

and 

X 3(iffp)t(l - f-2)-f - <f>c)/{zc - 2)K {A30) 
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in which z^" i has been assumed to be finite. 

Evaluation of the other Green's functions, Gi, G2 and G3, entails 
calculations which have many steps in common with that of Go- Since the 
integrals in (34), (42) and (52) differ from that in (16) only in that their 
integrands respectively contain the extra factors n, £uid n x e^, they can 
be rewritten as integrals of the form (A3) in which the functions 

fi(z/) = n/o, f2(z/) = ey,/o and £3(1/) = n x e^/o (^31) 

replace the S^{v) given by (A4). 

If po go ar^ correspondingly replaced, in accordance with (A5) and 
(A6), by 

P. - |(fi.|_e. +f^L=-c.)' fe- 1,2,3, (^32) 

and 

= |ci-i(ffc|^^^^ -f*L=-c.). ^ = 1'2,3, (^33) 

then every step of the analysis that led from (A7) to (A8) and (A9) would 
be equally applicable to the evaluation of Gfc. It follows, therefore, that 

Gfc'° ~ 2ci-^(l - [Pfc cos(| arcsinx) - Ciqfesin(| arcsinx)], Ixl < 1, 

(^34) 

and 

Gfc*""* ~ c-r^'bd - l)-»[pibsinh(iarccosh|x|) 

+ ciqfcsgn(x) sinh(|arccosh|x|)], Ixl > 1. {AZb) 

constitute the uniform asymptotic approximations to the functions G* inside 
and outside the laifurcation surface (the envelope) |x| = 1- 

Explicit expressions for and as functions of (r, z) may be fovmd 
from (8), (A19), and (A31)-(A33) jointly. The result is 

qi 

(A36) 

+ [kZ^^ ±R+^ + A^(^:^ =F^;^)]ev,i>}, (^37) 
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and 

-2-^(a;/c)(ffp)-^Cx=^5A-4{_(ip_i)[^:i ±^;t 

where use has been made of the fact that = — sin(^ — (pp)Brp + cos{(p — 
ipp)^tpp> Here, the expressions with the upper signs yield the and those 
with the lower signs the q^. 

The asymptotic value of each G^^"* is indeterminate on the bifurcation 
surface (the envelope). If we expand the numerator of (A35) in powers of its 
denominator and cancel out the common factor (x^ ~ 1)^ prior to evaluating 
the ratio in this equation, we obtain 

= ^^"''\^±i - (Pfc ± 2ciqO/(3ci2). (^39) 

This shows that Gfc^"*|^=^_ and Ga:*^"*U=^+ remain different even in the 
limit where the surfaces <j) == ^_ and <f> — coalesce. The coefficients qjk 
that specify the strengths of the discontinuities 

G'^°"%=^,-G*'^"*U=^_~t<l'tM (^40) 

reduce to 

qi ~ ^(a;/c)(ffp)-3[(l - |f2)fpe., + {Bp - z)e,^], (A41) 

q2~2i(a;/c)(ffp)-^e^^, (^42) 

and 

qa -2i (a;/c)(frp)-2[(zp - z)er^ - rpe^^] (A43) 

in the regime of validity of (A27) and (A28). 

When 0 < Zc - z «C (r^ — l)^rp, the expressions (A41) and (A43) 
further reduce to 

qi ~ ^(a;/c)(ffp)-2iii, and qa ~ 2t (a;/c)(ffp)-^n3, (^44) 

with 

m = (f-^ - ir)er^ - (1 - f~^)^e^P and ng = (1 - f-^)^erp + r-^e^^, 

{A45) 

for in this case (12b) — with the adopted plus sign — can be used to replace 
z-zp by (f^ - l)^fp. 
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